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We further study the previously proposed Ansatz, Tr(M^)=0, for a prediagonal hght Majorana 
type neutrino mass matrix. If CP violation is neglected this enables one to use the existing data on 
squared mass differences to estimate (up to a discrete ambiguity) the neutrino masses themselves. 
If it is assumed that only the conventional CP phase is present, the Ansatz enables us to estimate 
this phase in addition to all three masses. If it is assumed that only the two Majorana CP phases 
are present, the Ansatz enables us to present a one parameter family of solutions for the masses 
and phases. This enables us to obtain a simple "global" view of lepton number violation effects. 
Furthermore using an SO (10) motivation for the Ansatz suggests an amusing toy (clone) model in 
which the heavy neutrinos have the same mixing pattern and mass ratios as the light ones. In this 
case only their overall mass scale is not known (although it is constrained by the initial motivation). 
Using this toy model we make a rough estimate of the magnitude of the baryon to photon ratio 
induced by the leptogenesis mechanism. Solutions close to the CP conserving cases seem to be 
favored. 

PACS numbers: 



I. INTRODUCTION 

Remarkably, the recent KamLAND 0, SNO Q and K2K Q experiments have added so much to the results 
obtained from earlier solar neutrino, atmospheric Neutrino and a" elerator experiments i that our knowledge about 
the neutrino masses and presumed lepton mixing matrix is almost as great as our knowledge of the corresponding 
quantities in the quark sector. Still there is an uncertainty about the interpretation due to the results of the LSND 
experiment However, this experiment will be checked soon by the miniBoone collaboration so one can wait for 
confirmation before considering whether there is really a problem with the usual picture of three massive neutrinos. 
In any event, there is a strong presumption that this knowledge will play an important role in going beyond the 
standard model of electroweak interactions. 

One detail is, of course, lacking compared to the quark case. Since the neutrino oscillation experiments measure 
only the differences of the neutrino squared masses, the neutrino masses themselves are not known. According to the 
latest analysis j^j the best fit to these differences is: 

ml-mj = 6.9 x lO'^eV^, 
\'m\-m\\ = 2.6 X IQ-^eT^^. (1) 

Now, there is a simple complementary Ansatz for the 3x3 neutrino mass matrix, which, with some assumptions, 
enables one to obtain the neutrino masses themselves from Eq. (Q; it requires: 

Tr{M,) = 0. (2) 
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It should be remarked that is to be regarded as the prediagonal neutrino mass matrix. Furthermore, in the 
relation mi + m2 + ms = which evidently results if the neutrino mass matrix is taken to be real symmetric, the 
individual masses may be either positive or negative. The negative masses can be converted to positive ones by adding 
appropriate factors of i in the diagonalizing matrix. 

Eq. lI2Jl was motivated in from the grand unified model, SO(IO) Q and in by noting that it would hold if 
Ml, is the commutator of two other matrices, as may occur in certain models. If CP violation is neglected there are 
essentially two possible solutions of the Ansatz: either mi and m2 have the same sign and are approximately equal 
to each other and to — m3/2 or else TOi and TO2 have the opposite sign and are approximately equal to each other in 
magnitude but much larger than the magnitude of m^. 

In the present paper we will take the point of view that the Ansatz, Eq. is motivated from SO(IO). However, 
the analysis is of course not dependent on the motivation. The SO(IO) motivation arises from the observation that 
Eq. Q is, although it seems at first different, essentially the same as the characteristic prediction of grand unification: 

mi, = rmr, (3) 

relating the mass of the b quark with the mass of the the tau lepton (r « 3 takes account of the running of masses 
from the grand unification scale to the low energy hadronic scale of about 1 GeV Note that in SO(IO) the 

neutrino mass matrix takes the form: 

= Ml - mIm^HId, (4) 

where Ml, Mh and Md are respectively the mass matrices of the light neutrinos, heavy neutrinos and heavy-light 
mixing (or "Dirac matrix"). To start with, M^ is an arbitrary symmetric matrix. If it is real we have CP invariance. 
Generally the second, seesaw term is considered to dominate. However, as explained in 0, the present model 
is based on the assumption that the first term dominates. That might not be unreasonable since a rough order of 
magnitude estimate for the second term would be to^/IO^^ or about 3 x 10~^ eV. (The quantity 10"'^'^ includes a factor 

w 10). Thus the second term could be negligible if neutrino masses are appreciably greater than this value. 

In 7] the complementary Ansatz was mainly studied for the case of real M^. Here we will be primarily interested 
in the more general complex case which allows for non zero CP violation. Furthermore, the input squared mass 
differences were not taken to be very similar to those in Eq.(^ but were based on a least squares fit 12] of many 
different experiments including LSND. Here we shall adopt the more conventional values given in Eq.Q. A related 
analysis of Eq. was recently made in 1^ . 

For an understanding of the interesting leptogenesis mechanism ,^14] of baryogenesis it is important to also study 
the properties of the heavy neutrinos which appear there. In the present SO(IO) motivated framework this task turns 
out to be remarkably simple; the heavy neutrino mass matrix is given by 

Mh = cM,, (5) 

where c is a numerical constant. This means that the eigenvalues of Mh, to be denoted as Mi, M2, M^ are simple 
multiples of the light neutrino masses mi,m2,m3. In addition the unitary matrix, U which brings Mi, to diagonal 
form via 

U'^MiyU — diag{mi, 1712, ray,) — M^, (6) 

also diagonalizes Mh- In other words the heavy neutrinos are clones of the light neutrinos in this picture. The result 
follows from the choice of Higgs fields in SO(IO). Trilinear Yukawa terms which supply fermion masses can contain 
Higgs fields in the 10, 120 and 126 dimensional representations. To get the result just mentioned we need to require 
that there is only one "126" representation present, although any number of "10" 's and "120" 's are allowed. Of course 
we are also assuming the second term in Eq. to be negligible for the purpose of generating the light neutrino 
masses. 

In section II, we give our conventions for the lepton mixing matrix, including one conventional and two Majorana 
type CP violation phases. An approximate equation relating the complementary Ansatz to the parameters of the 
mixing matrix and the physical light neutrino masses is written down. The solutions for the neutrino masses in the 
CP conserving case, based on the results of neutrino oscillation experiments, are reviewed. In section HI, the Ansatz 
equation is solved on the assumption that only the conventional CP phase, 5 is non zero. It is found that the only 
solutions correspond to maximal phase, sin^6 = 1 and neutrino masses close to the ones obtained in the CP conserving 
case. In section IV we investigate the more complicated but very interesting case when only the two Majorana CP 
violation phases are non zero. In this case there is a family (modulo a discrete ambiguity) of solutions. We choose 
the mass of the third light neutrino, ma as our free parameter and calculate the remaining neutrino masses and the 
Majorana phases as functions of m^. The model gives a lower bound for and the cosmology criterion on the sum 
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of neutrino masses effectively yields an upper bound. The results for the full range are scanned numerically and a 
simple analytic interpretation of the pattern is presented. The neutrinoless double beta decay parameter, |mee| is also 
calculated for each value of m^. In section V we make a rough estimate of the baryon to photon ratio based on the 
leptogenesis mechanism. In order to do this it is necessary to make some statements on the masses and mixings of the 
heavy neutrinos. Our motivation for the original Ansatz suggests a "clone" model in which the heavy neutrinos have 
the same mass ratios and mixing matrix as the light ones. The only new parameter besides 7713 is the overall mass 
scale, which however is constrained by the original motivation to be somewhat on the large side. Then it is relatively 
easy to calculate the lepton asymmetry parameters, for the heavy neutrino decays as functions of mainly m^. We 
combine these quantities in a semi-quantitative way with criteria from previous treatments of the Boltzmann evolution 
equations for the decaying neutrinos. It is found that that the most plausible scenarios for leptogenesis involve small 
CP violating Majorana phases and light neutrino masses close to the ones predicted for the CP conserving cases.. 
Finally section VI contains a brief discussion and a brief summary. 



II. RELATING THE ANSATZ TO EXPERIMENT 



Here we will obtain an approximate equation which will be useful for relating the complementary ansatz to exper- 
imental information on neutrino squared mass differences and mixing angles in the general case where CP violation 
is allowed. The notation is the same as in section III of which contains more details. For convenience, we will use 
what seems to be the most common convention for the part of the leptonic mixing matrix, Kf.xp, which is measured in 
the usual neutrino oscillation experiments. This part can be constructed as a product of elementary transformations 
in the (12), (23) and (13) subspaces. For example in the (12) subspace one has: 



Wi2(.t:'12, 012 



) = 



COSf 12 

-e-''^i^sm6li2 




^sin9i2 







(7) 



with clear generalization to the (23) and (13) transformations. 
Then the usual convention corresponds to the choice: 

Kexp = W23(6'23, 0)^13(6*13, -(5)wi2 (6*12 , 0) , 

with three mixing angles and the CP violation phase S. Multiplying out yields: 



exp 



C12C13 
-S12C23 - Ci2Si3S23e"' 
S12S23 - Ci2Sl3C23e"^ 



S12C13 

C12C23 - Sl2Sl3S23e'' 



-C12S23 - Sl2Sl3C23e C13C23 



iS 



C13S23 



(8) 



(9) 



where Sij — sinOij and aj — cosOij . Since the neutrinos are of Majoran a type i n this model, there are expected 
also to be physical CP violating effects due to the Majorana phases [la llS llfl 13 • These may be introduced via a 
unimodular diagonal matrix of phases. 



T2+T3 = 0. 



The full lepton mixing matrix is then expressed as, 

K = KexpUJo^ir), 



(10) 



(11) 



which has three mixing angles and three independent CP violating phases. We shall use this form in what follows. As 
an aside, though, we remark that the full matrix could also be written P| in an unconventional, but more symmetrical, 
way as: 



K = a;23(^'23, 023)t^l3(6'l3, 013 )wi2 (6*12, </'l2)- 

As a final preliminary we need the leptonic W interaction term : 



(12) 



C = -^W;eLit.Kv + H.c., 



K 



(13) 
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where U is defined in Eq.(|SJl and fl^ is a unitary matrix which is needed to diagonahze the charged lepton mass matrix. 
At this point we shall make the common approximation that Q can be replaced by essentially the unit matrix. This 
is certainly not perfect but it seems reasonable for a start. Then U may be replaced by K , for which some elements 
are already well known. This enables us to present the ansatz in the form: 

Tr{M,K-lK:,j,iUo{2T)) = 0, (14) 

where Eqs. ((BJ and were used. With the parameterized mixing matrix of Eq. Q the ansatz reads: 

mie^^^i [l - 2i{ci2Si3fsinSe-'^] + 
TOae^*^^ [1 - 2i{si2Si3fsin6e''^] + 

mae^*^^ [i _,_ 2i{si3)'^ sinSe'^] = 0. (15) 

In this equation we can choose the diagonal masses mi,m2,7Ti3 to be real positive. However it will be a little more 
convenient in the CP conserving case to allow some of them to be negative as well as positive. We shall, for definiteness, 
mainly use the following best fit values for the mixing angles: 

s?2 = 0.30, s^3 = 0.50, s?3 = 0.003. (16) 

It should be remarked that the precise value of S13 is not well known, in contrast to the other two. 

Equation (|15|l contains three unknown masses and three unknown CP phases. It can be written as two real equations 
and augmented by two equations for two neutrino mass squared differences. Thus there are four equations for six 
unknowns. By assuming some special simplifications we can make the analysis tractable. 

For orientation let us first review the case when the theory is CP conserving so that all the three independent CP 
phases vanish. Then we will have 3 equations for 3 unknowns. The ansatz now reads mi + m2 + — 0. Define: 

A = (m2)2 - (mi)2 

B = (m3)2-(m2)2 (17) 

It can be deduced [w] from the experimental data that A is positive while the sign of B is not yet known. Their 
magnitudes are given in Eq. ^^1. Thus there are two separate cases to be considered. First consider both A and B 
positive. Then solving as in 7] gives the type I solution: 

mi = 0.0291 eV, to2 = 0.0302 eV, TO3 = -0.0593 eV. (18) 
Next consider the type II solution where B is negative; it gives: 

mi = 0.0503 eV, ma -0.0510 eV, m3 = 0.00068 eV. (19) 

Here mi and m2 are still almost degenerate but differ in sign. However m3 is now relatively small compared to the 
others. 



III. CONVENTIONAL CP VIOLATION 

A fully predictive simple case would correspond to keeping 6 as the only CP violation phase. Then the real part of 
the ansatz equation, (|15|l reads 



mi + m,2 + m3 — 2 sl^sin^S {0121^1 + 812^2 + m^) = 0, 
while the imaginary part yields 

Si3sm2i5(ci2mi + 812^2 — W3) = 0. 



(20) 
(21) 
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Note that the miS are being taken real here, ahhough they will be allowed to be either positive or negative. A 
negative is not a source of CP violation even though it corresponds to a Majorana phase of tt/2 when the masses 
are taken positive 20]. Now Eas. H17|l . H2()(l and 1)21(1 constitute 4 equations for the three to^'s and S. 

However it turns out that, except for the special case when Si^sin25 = 0, there is no consistent solution of this set 
of four equations for four unknowns. To see this, first consider solving simultaneously the three equations (|17|l and 
(|21|l when the special case does not hold. The numerical solution is seen to require B < and is found to be (with 
A>0): 

mi = -0.0548 eV, TO2 = 0.0554 eV, mg = -0.0217 eV. (22) 
We must now check to see if this is consistent with the remaining Eq. 1201) . That leads to the requirement: 

5^3^(5 =1(1+ "'^+"'" ) ^ 0.25, (23) 
4 TOa 

which, given the numerical value of s^g in Eq. (|16|l clearly leads to the contradiction siri^S « 80. This contradiction 
will persist even if the upper bound (about 0.044) rather than the best fit for sfg is used. The result is also not 
changed if the signs of all the mi's are reversed. 

Thus the only possibility for pure 5 type CP violation in the present scheme is the special case where sin^d = 1. 
Then we must solve simultaneously the three equations consisting of Eq. (|20() in which this substitution has been 
made for sin^S as well as Eqs. ()17|l. This results in the equation for, say TO2, 

stgn{mi)iml - Ay/^{1 ~ 2s%cl^) + TO2(1 - 2s%sl^) + [B + mlf'^l - 2.sl^) = 0, (24) 

where sign(m3) has been arbitrarily taken positive. Knowing m2, the other two masses may of course be obtained 
from Eqs. l(T7)l . 

Taking, for definiteness, the mixing angle from Eq. ((16(1 . one finds essentially two different solutions. These are 
quite similar to the Type I and type II solutions given above in the CP conserving case. The type I solution, with 
B > is 

mi = -0.0289 eV, m2 = -0.0301 eV, mg = 0.0592 eV. (25) 

The type II solution, with B <0, reads 

mi = 0.0503 eV, ma = -0.0510 eV, mg = 0.00081 eV. 

The very close similarity between the CP conserving solutions and the solutions with sin'^S 
due to the small value of s?q. 



IV. CP VIOLATION DUE TO MAJORANA TYPE PHASES 

Since, as we have just seen, there is only one particular allowed value for the conventional CP phase, S if it is 
considered as the only source of CP violation in the present scheme, it is of great interest to investigate the Majorana 
phases. Clearly, it seems sensible to study these phases with the simplification of putting 5 to zero. From Eqs. (O 
and it is seen that the same effect is accomplished by setting sig — 0. Then the ansatz equation l[T5]l takes the 
form 

mie^^^i + 77126^'^^ + mge^''^-' = 0. (27) 

For our present case it is convenient to take all three m^'s to be real and positive ( note that a phase angle = tt/2 
corresponds to what was taken as a shorthand to be a negative value of m^). Together, Eq. ((27(1 and Eq. 1(17(1 
comprise four real equations for five unknowns (three masses and two independent r^'s). To proceed we shall thus 
assume a value for mg so that we have four equations for four unknowns. In addition there is the two fold ambiguity 
due to the unknown sign of B. Finally we shall allow mg to vary to obtain a global picture of the situation. 

Now, once we have assumed a value for mg, we can immediately find mi and m2 from Eqs. ((17(1 . Furthermore, 
Eq. 1(27(1 can be pictured in the complex plane as a triangle formed from vectors with lengths m^, having angles 2ri 
as measured from the positive horizontal axis. Then let /i^ be the interior angle opposite side m^ as illustrated by the 
choice of triangle in Fig. ^ 



(26) 

= 1 is understandable 
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FIG. 1: Vector triangle representing Eq. 12711 . 



The problem reduces to one from elementary plane geometry. Given the three sides (m^), of a triangle, find the 
three interior angles (/i^). We may start for example, by using the law of cosines to get 

cos/ii = i- 2 3^ ^28) 

and continue similarly to get the others. Finally the parameters Ti which appear in the actual parameterization of 
Eq. (|10|) are found from Fig. 1 as 

n = g(7r - fii- 2^2), 
6 

= T(-27r-/ii +^2)- 
6 



(29) 



In particular, the quantities 



szn[2(ri - T2)] = -sm(/^i + ^12)^ 
sin[2{Ti - Ta)] = sm/Z2, 
sin[2{T2 - T^)] = -sinfii, 



(30) 



will turn out to be of interest. Actually, given the three interior angles /ii of a triangle we do not get a unique choice 
of phase differences (r^ — tj). While a rotation in the plane of the triangle will not change these phase differences, 
it is straightforward to see that the reflection of the triangle about any line in the plane will reverse the signs of all 
the phase differences. Thus there is another solution in which an extra minus sign appears on each right hand side of 

Now let us discuss the solutions of the complementary ansatz equation for various assumed values of TO3 . In Table 
^ the three real positive masses as well as the corresponding values of the two independent internal angles fii , /i2 of 
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the triangle are listed. Of course, /ii + /i2 + A'a = tt. The solution with B > (type I with ma > 1122) will be listed 
when it exists as well as the type II solution (i? < or TO2 > 1^3) ■ 

Let us start with large values of 7713 and go down. Just from the ansatz there is no upper bound on the value of 
7713. However there is a recent cosmology bound .21] which requires, 

+ |m2| + |to3| < 0.7 eV. (31) 

Thus values of ma greater than about 0.3 eV are physically disfavored. Tabled shows that at this value both type 
I and type II solutions exist and correspond to almost equilateral triangles. This is true also for higher values of 
7713. Notice that since the triangles are close to being equilateral, they have large interior angles and hence ( see for 
example Eq. (jSOl) large CP phases. The picture remains very similar down to around = 0.1 eV but as one gets 
closer to the value, roughly 0.0593 eV, where the real type I solution of Eq. H18|l exists, there is a marked change. 
It is seen that the interior angles of the type I solution become small as it prepares to go to the degenerate triangle 
corresponding to the real solution. We may get as small CP phases as we like by tuning close to the real solution; 
this is illustrated in Table HI for a particular value of 7713. If one further lowers m^, it is found that the type I solution 
no longer exists. On the other hand the type II solution persists and does not change much until 7713 approaches the 
small value of roughly 0.00068 eV. There are no solutions for 7713 smaller than this value. We can also tune 7713 as 
illustrated in the table to get as small CP phases as we like for the type II case. It should be remarked that the 
precise numbers in Tableware based on the assumption that the best fit numbers given in Eq. are exact and hence 
are meaningful to the accuracy given only in the sense of comparing the various solutions with each other, not with 
experiment. It is straightforward to give an analytic interpretation of the pattern of solutions just observed. First 



type 7771,7712,7713 in eV MiiMa in radians |777ee| in eV ei, £2, £3 



I 





2955, 





2956, 





3000 


1 


038, 1.039 





185 





342, 0.433, 0.017 


II 





3042, 





3043, 





3000 


1 


055,1.056 





187 





330, 0.426, -0.0172 


I 





0856, 





0860, 





1000 





946, 0.952 





058 





138, 0.060, 0.00137 


II 





1119, 





1122, 





1000 


1 


106, 1.111 





065 





194, 0.088-0.0024 


I 





0305, 





0316, 





0600 





258, 0.268 





030 





00982, 0.00422, 0.00004 


II 





0783, 





0787, 





0600 


1 


172, 1.186 





043 





094, 0.041,-0.0011 


I 





0291, 





0302, 





0592715649 





000552, 0.000574 





030 


1 


96 xl0~^ 0.84 xl0~■^ 0.71 xlO"'^ 


II 





0774, 





0782, 





0592715649 


1 


174, 1.188 





042 





047, 0.020, -0.0011 


II 





0643, 





0648, 





0400 


1 


243, 1.268 





033 





052, 0.023,-0.000681 


II 





0541, 





0548, 





0200 


1 


355. 1.417 





024 





018, 0.0078,-0.000335 


II 





0506, 





0512, 





0050 


1 


386, 1.658 





021 





0057, 0.0025,-0.0000824 


II 





0503, 





0510, 





0010 





814, 2.313 





021 





00073, 0.00031,-0.0000122 


II 





0503, 





0510, 





0006811 





051361, 3.089536 





021 





0000348, 0.0000150, -0.601 xlO"® 



TABLE I: Panorama of solutions as ma is lowered from about the highest value which is experimentally reasonable to the 
lowest value imposed by the model. In the type I solutions 7713 is the largest mass while in the type II solutions 7773 is the 
smallest mass. For each value of 7773, the values of the model predictions for 7771 and 7772 as well as the internal angles ^1 
and jj,2 are given. The model prediction for the neutrinoless double beta decay quantity |mee| is next shown. Finally, the last 
column shows the estimated lepton asymmetries due to the decays of the heavy neutrinos. Note that the reversed sign of lepton 
asymmetry is also possible, as discussed in the text. 

note that CP violation corresponds to a non degenerate triangle. Note also that the orientation of the triangle in the 
complex plane is just obtained by imposing the unimodularity condition for w(t) in Eq. H10|l . Hence the internal 
angles, are really the intrinsic carriers of CP violation. The determinant of whether one has CP violation is the 
non vanishing of the quantity: 

A — [777(777 — 777l)(777 — 7772)(777 — 7773)]"'"''^, 

771 = -(7771+7772 + 7773), (32) 

which just expresses the area, ^ of a triangle in terms of the lengths of its sides. This area may be rewritten in the 
convenient form: 

A^J ([(7771 + 7772)2 - 777^] [777^ - (777l - 7772)^]) . (33) 

Now we may see that the vanishing of the first factor corresponds to the type I real solution while the vanishing of 
the second factor corresponds to the type II real solution. Furthermore, for a solution to exist, the argument of the 
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square root should be positive. With the second factor, that estabhshes the minimum ahowed value of while with 
the first factor, that establishes the minimum value of which allows a type I solution. 

An important test of the model is the experimental bound on neutrinoless double beta decay. This implies 

|mee| < (0.35 ^ 1.30) eV, (34) 

where 

\mee\ = \mi{Ke.piife-^'^' +m2iKercpl2fe-^'^' +m3{Ke.pl3fe-^'^'\. (35) 

Using the parameterization of Eq. (jSJ and approximating S13 = (which is reasonable in the present model since 
is never much larger than mi or m2), this can be written simply as: 

|mee| = [(Ci2"ll)^ + (s?2"t2)^ + 2mim2Ci2Si2COs(/ii + A*2)]"^^^- (36) 

Here Eqs. (|29|) were also used. Reading from then enables us to calculate |mee| for each line of Table I. It is 
seen that \mee\ decreases smoothly with decreasing 7713 for each of the type I and type II solutions. All the values of 
ma listed are consistent with the present bound. It is interesting that an improvement of the experimental bound by 
an order of magnitude j'SsI would provide a good test of the model. 



V. ESTIMATE FOR LEPTOGENESIS 



When one adopts the SO(IO) motivation for the present Ansatz, it turns out that the resulting model predicts in a 
simple way the properties of the heavy neutrinos which are intrinsically contained in the SO(IO) theory. This feature 
may be used in connection with the leptogenesis mechanism |14| of baryogenesis. According to this mechanism, the 
CP violating and lepton number violating decays of the heavy neutrinos at a high temperature (corresponding to the 
grand unification scale) in the very early universe establish a lepton asymmetry. As the universe cools further, the 
(B+L) violating but (B-L) conserving "sphaleron" interaction j2J] converts this into a baryon asymmetry which may 
be compared with the observed ratio of baryons to photons in the universe. There are many interesting discussions of 
this mechanism in the literature |25l l27l |28L l29l l30l |. Here, we will estimate the dependence on neutrino masses and 
CP phases of the predicted baryon asymmetry in the present model. 

The starting point of this discussion is the Yukawa term of the Lagrangian density which describes the tree level 
decay of a heavy Majorana neutrino, Nj (where the subscript denotes a 3- valued generation index) to a Higgs doublet 
member, 



plus the appropriate member of the left-handed lepton doublet. 



Then the Yukawa term reads: 



L.-i'^]- (38) 



Cvukawa - - ^ L,X,,<i>^N, + H.C., (39) 



where A^j is the matrix of Yukawa coupling constants. We can simplify this expression, which is supposed to contain 
the fermion fields in prediagonal bases, in several ways. First, at the high temperature for which the N decays are 
relevant, the phase transition to spontaneously broken SU(2) x U(l) has not yet taken place. Thus we can consider 
the light fermions in Li to be massless and there is no need to insert suitable unitary matrices to bring the light field 
mass matrices to diagonal form. However the heavy neutrino, N should be related to the physical field N with a 
unitary matrix U a.s N = UN. As mentioned in section I, if the SO(IO) model contains only a single "126" Higgs 
type field (although any number of "10" 's and "120" 's are allowed) and also if the first (non see-saw) term in Eq. © 
is dominant, the prediagonal mass matrices for the light and heavy neutrinos must be proportional to each other and 
the diagonalizing matrix U must be the same one which appears in Eq. (|13() . Approximating, as we did earlier, fl to 
be essentially the unit matrix we can set U K. li the model of section IV is adopted, for example, we can specify 
K, including CP phases, to a fair approximation for each assumed value of 7713. Finally we approximate the matrix 
of Yukawa couplings by: 
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where mf are the three charge 2/3 quark masses at a low energy scale, r' « 3 is a suitable factor for running these 
masses from the grand unified scale to the low energy scale and < 0*^ >w 246/ \/2 GeV. Note that Cvukawa is the 
term responsible for generating the neutrino Dirac matrix, Md in Eq. In the simplest approximation to the 

SO(IO) theory the charge 2/3 quark mass matrix and neutrino Dirac mass matrix are proportional to each other and 
diagonal (since the quark mixings are after all small). 

Putting these things together we arrive at the "effective" term for calculating the heavy neutrino decays (at grand 
unified scale temperature): 



^Yuka 

where 



^-J2L^h,,<^>^N,+H.c., (41) 
■' < 4P > r' 

The quantities needed for the calculation are the matrix products {h'^h)ij. We may further simplify these products 
by noting that the top quark mass is much heavier than the others so the products approximately become hl^h^j. 
Specifically, the diagonal products are: 

{h^h)n - {si2S23/r')\ 

ih^h)22 ~ (Cl2S23/r')^ 

(h^hhs « (cl2/r')^ (43) 

where we used the numerical coincidence that mt —< (j) >• Furthermore, we set S13 = in agreement with the 
model of section IV ( See the parameterization of Eq. @). Numerically, with Eq. (|16() and (r')^ « 10 one obtains 
{K^h)ii w 1.50 X 10"^ {h^h)22 ~ 3.50 x lO'^ and {h^h):^^ w 7.00 x 10^^^ jj^ ^gj-^^s of these diagonal products, the 
tree level widths of the heavy neutrinos are given by, 

r, = <il^. (44, 

where Mi is the mass of the ith heavy neutrino. The off diagonal products play an important role in determining the 
lepton asymmetry. They are explicitly given in the model of section IV as: 

{h^h)i2 ~ -si2Ci2sl-,e'^^'-^-^/ir')\ 

{h^h)i3 « Sl2S23C23Cl3e^^"^""^V(^')', 
{h^h)23 ~ -S23Cl2C23Cl3e'^"^^"^V(^')', 

(h^h),, = (h^h)*,, (45) 

where the CP phases depend on the choice of TO3 as explained in section IV. Numerically, one has {h''h)i2 ~ 
-2.29 X 10-^exp[iiTi - T2)], (/it/i)i3 w 2.74 x 10~^exp[i{Ti - T3)] and (/it/i)23 « -4.18 x 10-^exp[iiT2 - T3)]. In 
arriving at these estimates from Eq. H16|l we arbitrarily took all the signs of the trigonometric functions to be 
positive. This will not lead to any ambiguity since, for the application of interest, the off-diagonal products must be 
squared. 

The lepton asymmetry e^, due to the decay of the ith heavy neutrino is defined as the ratio of decay widths: 

^ r(jv,^L + $)-r(jv.^Z + $) 

' r(iVj ^ L + $) + T{Ni ^ i + $) ■ ^ ' 

In this formula i + (f> stands for all pairs of the types ej + 0+ and + (fP . This is an effect which violates C and 
CP conservation, in agreement with the requirement of Sakharov |25| | . To get a non-zero value one must include the 
interference between the tree diagram from Eq. H41|l and the one loop diagrams (of both "self-energy" and "triangle" 
types). If the masses of the heavy neutrinos are well separated the result J^^^] is: 



Im[{h''h),j{h''h)iji /•/7i^2/,^2n 



+ 1 - (1 + j;);n(l -I- 1/x) 



(47) 
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Note that the contribution to the lepton asymmetry of the hghtest heavy neutrino is expected to be the most important 
one for the final calculation of baryon asymmetry |27j . 

Now let us make numerical estimates for the lepton asymmetries when all CP violation is due to Majorana phases 
(section IV). From Eq. I^l we relate the heavy neutrino masses to the light neutrino masses simply as: 

M, = cnii, (48) 

where c is a real, positive constant. This equation has earlier been used [s^ for the study of leptogenesis in the 
framework of a left-right symmetric model. It should be noted that renormalization group effects |33 will modify the 
exact proportionality of the light and heavy neutrino masses as well as the equality of the corresponding diagonalizing 
matrices. This should be taken into account for a more accurate treatment. In the model of section IV the third 
neutrino is typically somewhat further away in mass from the other two, which are always relatively close. For 
example, in the type II situation, 7713 is the lightest of the light neutrino masses so Afa will be the lightest of the heavy 
neutrino masses and the contribution to the lepton asymmetry is 63. Using Eqs. (|47|l . H48(l . 145|l and H30|l we obtain: 

£3 « [~A.27f{[mi/m3f)sinfi2 + 9.94/([m2/m3]^)sm^i] x 10~^ (49) 

Notice that c has canceled out in this formula and also cancels out in the determination of the angles /i^. Thus the 
lepton asymmetry given by this formula does not depend on the overall scale of the heavy neutrino masses. 

In the type I case, the heavy neutrino spectrum consists of two nearly degenerate lighter states, (A^i, N2) and a 
heavier state, N^. For the corresponding asymmetries ei and e2, the diagrams involving self energy type corrections 
are enhanced since an internal heavy neutrino line will be close to its mass shell. The formulas [^H thus, for greater 
accuracy, involve the decay widths and we will approximate: 

Im[{h^h)i2ih^h)i2] {Ml - Mi)MiT2 

ei - 



{h^h)ii{h^h)22 {Ml - M|)2 + MlTl ' 
_ Im[{h^h)i2{h^h)i2] {Ml - MDMaFi 

{h^h)i^{h^h)22 {Ml - Ml f + MlTl ^^^> 

Again, we may replace the heavy neutrino masses using Eq. (j48|l and note that the factor c cancels out. Inserting 
numbers, we obtain: 



ei 



1.39 X 10 '^sin{pi + ii2)'niim2{m\ — m\) 
{ml — m\Y + 1.94 x \{)^^m\m2 



''2 

5.96 X 10^^sm(/ii + /i2)"^i'Ti2(wi — ^2) 
^ {ml - jnlY + 3.56 x lQ-''m\ml ' ^^^^ 

The values of all these asymmetries for the range of possibilities are listed in the last column of table I. 

Furthermore it must be noted that, owing to the non-uniqueness of sign for all of Eqs. (|30|l . reversing the signs of 
all the lepton asymmetries also yields a solution corresponding to our initial Ansatz. 

Although the scale of the heavy neutrinos has been seen to cancel out of the formulas H49f) and H51|l for the lepton 
asymmetries in favor of their ratios (which are the same as those of the light neutrinos in this model), there is 
nevertheless a consistency condition implied by the SO(IO) motivation for the starting Ansatz. This arises since Eq. 
(|39|l is not only the source of the lepton asymmetry but also provides the seesaw contribution to the light neutrino 
masses. For our motivation we assumed that this contribution was dominated by the first term of Eq. I^J. To make 
a rough estimate of what this means we assume all matrices of the seesaw term to be diagonal. Then the value of 
c^/^ defined in Eq. (|48|l should be greater than mf / {m,ir') in order that the first term of Eq. Q be greater than the 
second term. In the case of the type I solution with m3 « 0.06 eV shown in Table I, this implies that the lightest 
heavy neutrino should be heavier than about 2.6 xlO^'^ GeV. For the case of the type II solution with m^ « 7 x 10"'' 
eV, the lightest heavy neutrino should be heavier than about 4.4 x lO'^^ GeV. 

The goal of the baryogenesis problem is to understand the ratio ris = ns/n^, the net baryon number density 
divided by the photon density. Experimentally, this quantity is found 21], from the study of big-bang nucleosynthesis, 
to be 

r^B = (6.5 ±0.4) X 10^1°. (52) 

To obtain non-zero rjB, it is not sufficient, as pointed out by Sakharov|25j. just to have non-zero values of the lepton 
asymmetry e^, defined in Eq. (|46|) . In addition, the CP violating decays of the heavy neutrinos must occur out of 
thermal equilibrium. A detailed treatment requires solution of the Boltzmann evolution equations for the system |34| . 
Here we shall make a rough estimate which we use to draw what might be a fairly robust conclusion. 
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First, we should remark that the baryon asymmetry generated by the sphaleron mechanism would be about -1/3 
|26j |. (for a review see 113) initial lepton asymmetry. The lepton number violating decays of the ith heavy 

neutrino are usually roughly taken to be out of equilibrium if the decay rate Ti in Eq. 1)44(1 is less than the Hubble 
rate, 

H w 1.7gl^^T^/Mp, (53) 

where « 100 is the number of effective light degrees of freedom at the leptogenesis scale, T is the temperature ( 
corresponding to the mass of the decaying heavy neutrino) and Mp 1.22 x 10^^ GeV. In the present model this 
ratio takes the explicit form: 

H 427Mi ' ^ ' 

which is seen to be inversely proportional to Af^. The net baryon asymmetry is estimated as j27j . 

VB ^ -:7—^e,D^, (55) 

where the Di are suppression factors to be obtained by numerical solution js^l of the Boltzmann equations. It is 
generally accepted that only the contributions of the lightest heavy neutrinos should not get washed out; thus we will 
set Di =0 for the heavier neutrinos. If 10 < Ki < 10^, the suppression factor is often approximated by the analytic 
form I35I 



0.3 



r.uJr..,0., - (56) 



When Ki < 1, the suppression factor is expected to be of order unity if Mi is not too large. However, as Mi gets 
larger there is a sizeable washout effect 36J . 

Glancing at the last column in table I and comparing with the experimental value of 775 in Eq. H52() as well as Eqs. 
((55|l and ((56|l suggests that the values of jys obtained for typical values of the assumed light neutrino mass parameter 
TO3 would be considerably larger than the experimental baryon asymmetry. However, we can expect to be able to 
obtain agreement with the experimental value since, as discussed in section IV, we may make the Majorana CP phases 
as small as we like by continuously tuning the independently chosen variable, 7713 so that the triangle of mass vectors 
gets arbitrarily close to one of the two degenerate straight line cases which causes A in Eq. to vanish. Thus the 
solutions of the model which would be consistent with the observed baryon asymmetry correspond to neutrino masses 
more or less close the real cases of either Eq. (|18|) or Eq. (I19|l . 

The qualitative points: i. that in the present model the value of the free parameter, 7713 can always be tuned to 
be arbitrarily close to its values for the two real solutions (so that the CP violation and hence leptogenesis strength 
becomes as small as desired) and ii. that the characteristic lepton asymmetries, for values of away from these two 
real solutions are rather large, comprise the main result of our discussion of the application of the Tr{Mv) ~ Ansatz 
to the baryogenesis problem. These points lead to the expectation that the physical value of is likely to be close 
to one of the two values in Eq. (|18|) or Eq. (|19|l and that this conclusion might persist even when our simplifications 
are not made. A more accurate treatment would include the features: a) effect of non-trivial charged lepton mixing 
matrix, b) renormalization group induced deviations from the "clone" treatment of the heavy Majorana neutrinos 
and c) full integration of the Boltzmann evolution equations. Even though our main conclusion is a qualitative one, it 
seems nevertheless an interesting exercise to find what values of light and heavy neutrino masses corrrespond to the 
correct order of magnitude of the observed baryon asymmetry. Note that the seeming great accuracy of the entries in 
Table I is not meant for precise comparison with experiment, but for comparison of the results of different 7TI3 choices 
with each other. 

Specifically, consider the tuned type I solution in table I with « 0.05927 eV. We noted in the discussion after Eq. 
that this would correspond to heavy neutrino "clone" masses (Mi, M2, M3) greater than about (2.60, 2.70, 5.27) x 
10^^ GeV, respectively. We assume that the two lighter neutrinos are the important ones and set D-^ = 0. The ratios 
{Ki^K2) defined in Eq. H54|) would then be less than about (16.5, 37.1) and would result in suppression factors 
{Di,D2) greater than (0.010, 0.0037). Using Eq. |(SSJ) and table I for ei and e2 then gives |77b| = 5.4 xlO~^°, close 
to the experimental value in Ea. H52ll . This can be adjusted by further tuning 7713 or to some extent by varying the 
overall mass scale of (Mi, M2, M3). 

For the type II case, first consider the solution in table I with 771,3 ~ 0.0007 eV. As discussed before, this would 
correspond to heavy neutrino "clone" masses (Mi, M2, M3) greater than about (320, 320, 4.4) xlO^^ GeV. In this case, 
M3 is the lightest of the three heavy neutrinos and is assumed to be the relevant one. We thus set Di — D2 = 0. The 
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ratio K3 given in Eq. (|54|l is then about 0.45 and indicates that the lightest heavy neutrino is, as desired, decaying 
out of equihbium. However, because its mass is considerably higher than that of the type I case just discussed, there 
is more wash out 33|, « 3.5 x 10~^. Reading £3 from table I then gives 775 ~ 5 x 10~^^, about three orders of 
magnitude too small. Thus we must raise the value of a bit. Backing off a little to the case 7713 « 0.005 eV in 
table I increases the value of €3 and also allows us (in line with the dominance of the first term in Eq.Q)) to choose 
the lower bound of M3 to be smaller, around 6 xlO^'' GeV. This results in an estimate \t]b \ ~ 16 x 10^^*^, which is 
the correct order of magnitude. One might wonder whether the contributions to 775 from ei and e2 are completely 
washed out in a case like the present. However, even if they were dominant, it would just require us to tune more 
closely toward small m^. 

Thus, if the model of CP violation with just the Majorana phases is correct, the magnitude of the baryon to photon 
ratio can be understood when either the sum of the three light neutrino masses is about 0.118 eV and \mee\ =0.030 
eV (type I) or the sum of the three light neutrino masses is about 0.107 eV and |7Tiee|= 0.021 eV (type H). In both 
cases the CP violating Majorana phases are extremely small. That might suggest a possible model in which a small 
CP violating perturbation due to some separate effect modifies an otherwise CP conserving lepton sector. 

We can also calculate the baryon to photon ratio in the model of section HI, where S is the only CP violating 
phase. There we noted that the only possible choices of S consistent with our Ansatz satisfy siri^S — 1. Then we 
have the type I solution for light neutrino masses given in Eq. (|25() and the type II solution given in Eq. (|26() . The 
corresponding CP violation factors are now (to first order in the small parameter S13) for the type I case: 



,n4 C23S23, 



(r') 

Im[{h^h)i3{h^h)i3] « S23C^3, 

Im[{h^h)23{h^h)23] - -Im[{h^h)i3{h^h)i3], (57) 



and for the type II case: 



T \thiu\ fuU\ 1 ^ si3sinSsin{20i2) 3 
Im[{h^h)i2{h^h) 12] w C23S23' 

T ffhth\ 1 ^ si3sinSsin{29i2) 3 

Im[{h^h)i3{h^h)i3] « S23C^3, 

Im[{h^h)23{h^h)23] - Im[{h^h)i3{h^h)i3], (58) 

As in the cases where only the Majorana phases contribute to the CP violation, the predicted lepton asymmetries, 
ei will typically lead to a value of the baryon to photon ratio much larger than the experimental one. In the present 
case it is not possible to fine tune S. The only possibility would be to fine tune S13 to an extremely small value. This 
seems more artificial since S13 is not required to vanish in the CP conserving situation. In any event, an experimental 
measurement of non-zero S13 would, practically speaking, rule out this case as a candidate for leptogenesis. 



VI. DISCUSSION AND SUMMARY 



In this paper, we investigated an Ansatz which correlates information about the four quantities in the light neu- 
trino sector which are not yet known from experiment; namely, the absolute mass of any particular neutrino, the 
"conventional" CP violation phase and the two Majorana phases. Of course, with input from analyses of neutrino 
oscillation experiments, the masses of the other two neutrinos can be found, up to a discrete ambiguity, if the mass of 
one is specified. The results of the present paper can be used for calculating many quantities of experimental interest 
like the neutrinoless double beta decay amplitude factor rriee (presented in section IV) and various lepton number 
violating decays. 

The Ansatz is not completely predictive, unless some assumptions are made. We first reviewed the case of assumed 
CP conservation (where just the three neutrino masses are obtained). Then we showed that if only the "conventional" 
CP violation phase is assumed to be non-zero, its value is fixed by the Ansatz to be maximal. A possibly more 
interesting case appears if we assume that only the two Majorana phases are non-zero. This enables us to scan the 
limited allowed range of assumed neutrino mass, 7713 (say) and find the other two neutrino masses as well as the two 
Majorana phases for each value of 1713. The result seems to cut through a "cross section" of interesting possibilities 
which are described in a simple way. The still more complicated case without setting any of the three CP phases to 
zero gives a two parameter family of solutions and will be treated elsewhere. Another (common) assumption we made 
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for a first analysis is that the measured lepton mixing matrix is dominated by the neutrino factor. This is consistent 
with the finding in recent years that the mixing in the neutrino sector is apparently much larger than the mixing in 
the quark sector (which in models is usually relatively small and similar to that of the charged lepton factor). 

It seems relevant to discuss briefly the status of the motivations for the complementary Ansatz we are using. One 
motivation, based on a loop mechanism for generating neutrino masses was discussed recently by He and Zee j^. 
Our motivation 1^ was based on the grand unification group SO(IO). This group is well known to have the elegant 
feature that it accomodates one generation of elementary fermions as well as an extra (now desired) neutrino field 
in its fundamental spinor irreducible representation. Naturally, there are many possibilities for doing a detailed 
calculation using this group. One may ask whether it should be regarded as being derived from a superstring theory, 
whether it should be supersymmetric, whether the symmetry breakdown should be dynamical, whether the symmetry 
breakdown should be induced by Higgs fields and if so what kind and how many, etc? We focus, in our motivation, 
on the conventional possibility of using Higgs fields since it seems almost kinematical now (although since no Higgs 
field has yet been seen one should keep an open mind). Of course, there have been many interesting treatments along 
these lines [2,7^ . Our Ansatz is suggested by a relation involving only the neutrino mass matrix which might be true 
(or at least approximately correct) in a large number of models. In SO(IO), tree level masses from a renormalizable 
Lagrangian can be obtained by using any number of Higgs mesons belonging to the 10, 120 and 126 dimensional 
representations. However, examining the form of the predicted mass matrices shows that the following fairly general 
relation ^ holds: 

Tr{M° - rM^) cx Tr{ML), (59) 

for any number of lO's and 120's but only a single 126 present. Here and are respectively the prediagonal 
mass matrices of the charge -1/3 quarks and charge -1 leptons while r w 3, as previously mentioned. M^, which arises 
from the 126 Higgs field Yukawa couplings, is the non seesaw part of the light neutrino mass matrix which appears 
in Eq. Taking traces cancels the contributions (antisymmetric matrices) of any 120 Higgs multiplets to the left 
hand side. Then, assuming the transformations which bring and to diagonal form to be roughly close to the 
identity we observe that the left hand side is approximately equal to nib — rnir, which is about zero. In fact this is 
a characteristic prediction of grand unification. In turn the right hand side gives us the starting Ansatz when it is 
assumed that the non seesaw term dominates in Eq. Of course, if this domination is to hold the masses of the 
heavy neutrinos should not be too low. The present paper is in effect exploring the range of possibilities which exist 
when these assumptions are made in SO(IO) models. An interesting question is whether this kind of limit or the pure 
seesaw limit gives a better description of nature, even if both terms are actually required. 
We remark that SO(IO) also gives another similar relation, 

Tr{M^ - r'Mo) oc Tr{ML), (60) 

when only one 126 Higgs field exists. Here and AId respectively denote the prediagonal charge 2/3 quark mass 
matrix and the prediagonal neutrino "Dirac" matrix connecting the heavy and light neutrino fields. 

An intriguing way to learn more about CP violation in the lepton sector is the study of the leptogenesis mechanism 
of baryogenesis. We saw that the treatment of this process simplifies when one adopts the present SO(IO) motivation. 
Then the light neutrino mass matrix and are approximately equal and proportional (due to the assumption 
of only one 126 field in the theory) to the heavy neutrino mass matrix Mh- The only free parameter for the heavy 
neutrinos is their overall mass scale and this should not be too small to preserve non seesaw dominance. We showed in 
section V that it is easy to estimate the lepton asymmetry parameters ti for a "panorama" of values of the independent 
variable since they are actually independent of the overall heavy neutrino mass scale. As far as the resulting baryon 
to photon ratio, rjB (parameterized in Eq. \bb\ ) is concerned, the typical values of the e,; give rjB much greater than 
the experimental one for suppression factors Di of order unity. We observed that if the suppression factors are not too 
small one can therefore always choose a value of 7713 close enough to one of the two essentially different CP conserving 
solutions so that the Majorana phases are small enough to get experimental agreement for rjB- Using estimates of the 
suppression factors taken from other earlier studies, we noted that this conclusion seems reasonable. Of course the 
study of the suppression factors by solving the Boltzmann evolution equations is an important topic which involves 
many subtleties and would repay further work in the present model. Finally, the posible indication of very small CP 
phases might suggest a model in which the CP violation in the lepton sector has a separate identifiable source. 
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